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Abstract .  Theoretical  calculations  of  wall  pressures  and 
elastic  wave  radiations  from  high-explosive  detonations  in  spherical 
cavities  in  salt  are  compared  with  close-in  experimental  results 
from  the  Cowboy  Program.  Spherical  charges  of  Pelletol  ranging 
from  20  to  2000  pounds  (9  to  900  kg)  were  detonated  in  spherical 
cavities  12  and  JO  feet  in  diameter  (3.66  and  9-1^  meters).  The 
elastic  wave  and  the  wall  pressure  results  were  found  to  be  sensitive 
to  the  explosive  detonation  details.  New  ways  are  suggested  for 
correlating  such  quantities  as  explosive  and  radiated  energies,  peak 
values  of  particle  displacement,  velocity,  acceleration,  and  radial 
stress.  This  correlation  procedure  can  be  useful  in  scaling  peak 
particle  velocities  from  cavity  or  uncoupled  explosions,  and  is 
demonstrated  for  the  Cowboy  results. 


*A'ny  Vievs  expressed  in  this  paper  are  those  of  the  authors . 
They  should  ;\ot  be  Interpreted  as  reflecting  the  views  of  The  RATIO 
Corporation  c.r  the  official  opinion  or  policy  of  any  of  its  govern¬ 
mental  or  private  research  sponsors.  Papers  are  reproduced  by  The 
RAND  Corporation  as  a  courtesy  to  members  of  its  staff. 
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INTRODUCTION 


This  paper  reports  the  results  of  an  investigation  [Erode,  I96O; 
Parkin,  I962]  which  was  stimulated  by  the  theoretical  findin-gs  of 
A.  L.  Latter  and  his  collaborators  at  RAND  concerning  the  concealment 
of  underground  nuclear  explosions  [Latter  et  al.,  I959j  I961].  A 
central  idea  behind  this  method  of  concealment,  or  "decoupling,"  is 
that  if  a  detonation  occurs  in  an  underground  cavity  which  is  large 
enough  to  cause  all  displacements  in  the  surrounding  rock  to  stay 
within  the  "elastic"  range,  most  of  the  explosive  energy  should  leave 
the  cavity  in  the  form  of  seismic  waves  of  high  frequency  which  would 
not  radiate  to  large  distances.  A  very  considerable  fraction  of  the 
explosion  energy  should  also  remain  within  the  hot  gases  contained 
by  the  cavity  after  the  explosion.  Clearly,  the  difficulty  of 
detecting  and  identifying  such  a  "decoupled"  nuclear  explosion  would 
thereby  be  increased. 

These  ideas  have  since  received  some  experimental  verification 
from  the  underground  tests  of  project  Cowboy  [Herbst  et  al.,  19^1; 
Adams  and  Allen,  I961].  In  these  tests,  extensive  seismic  data  were 
obtained  when  high  explosives  were  detonated  at  the  center  of  two 
spherical  cavities,  of  6  ft  (I.89  meters)  and  15  ft  (^.57  meters) 
radii,  in  a  salt  dome  and  also  when  tamped  charges  were  exploded 
in  the  salt.  Comparison  of  seismic  signals  from  cavity  and  tamped 
explosions  showed  significant  decoupling-  In  addition  to  the  seismic 
data,  a  considerable  body  of  measurement  was  obtained  from  the  ca.vity 
explosions  which  record  the  close-in  effects  of  wave  propagation  in 
the  salt  ;ina  the  explosive  blast  pressures  on  the  cavity  wall 
^Mirphey,  I961]. 


Although  some  theoretical  calculations  of  events  in  the  salt 
near  the  cavities  have  been  compared  with  experiment  [Herbst  et  al, 
1961]  and  although  Independent  calculations  of  the  hydrodynamics 
inside  the  cavities  have  been  carried  out  at  Lawrence  Radiation 
Laboratory  [H-  R-  Nuckolls,  private  communication],  as  far  as  the 
authors  are  aware,  no  results  of  an  entirely  theoretical  nature 
concerning  both  cavity  hydrodynamics  and  close-in  phenomenology  in 
the  salt  have  been  published.  Starting  from  first  principles, 
theoretical  calculations  of  the  cavity  hydrodynamics  are  used  as 
the  basis  for  a  series  of  elastic  wave  calculations  in  order  to 
obtain  theoretical  descriptions  of  the  close-in  elastic  radiation 
field  in  the  salt  and  the  nonsteady  flows  in  the  cavities.  These 
theoretical  results  are  compared  with  experiment. 

The  nature  of  these  results  is  such  as  to  indicate  some  physical 
trends  which  might  otherwise  be  obscured  in  empirical  or  purely 
experimental  determinations.  These  calculations,  when  compared  with 
experiment,  permit  the  assessment  of  the  precision  of  both  elastic 
wave  theory  and  of  the  theoretical  hydrodynamic  description  of  the 
cavity  blast.  We  point  also  to  some  methods  of  correlating  the 
experimental  results  which  may  previously  have  escaped  notice. 
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CAVITY  EXPLOSION  THEORY 


Explosive  properties- -loading  density  and  equation  of  state 
effects .  The  initial  density  (loading  density)  of  the  explosive 
plays  an  important  role  in  determining  both  the  detonation 
properties  and  the  final  gas  properties.  Jones  ana  Miller  [l9'+8] 
give  compositions  and  equations  of  state  for  the  gaseous  products 
from  TNT  at  loading  densities  of  1.5  and  1.0  g/em^.  The  gaseous 
products  are  remarkably  different  both  at  the  Chapman- Jouget  point 
(at  the  detonation  front)  and  after  expansion  to  low  pressure. 
Pelletol  was  used  in  the  Cowboy  experiments.  It  consists  of  small 
spheres  of  solid  TNT;  which  results  in  an  over-all  loaning  density 
of  unity  (p^  =  1.0  g/cm^). 

The  detonation  wave  in  center-detonated  cast  TI’JT  ( =  1.5  g/cm^ 
leads  to  a  peeik  ( Chapman- Jouget )  pressure  of  about  l6o  kb,  an  energy 
density  of  about  5*5  x  10^^  ergs/g;  and  a  detonation  velocity  of 
6.4  m/msec.  This  results  in  a  temperature  close  to  and  a 

compressed  density  of  about  2.0  g/cm^. 

In  contrast;  the  similar  detonation  front  conditions  in  TNT  of 
loading  density  unity  (o^  =  l-O  g/c's?)  (Pelletol)  give  pressure 
of  about  82  kb;  an  energy  density  of  about  4  x  10^*^  args/g;  -i 
detonation  velocity  of  about  4.5  m/mseC;  a  temperature  of  ^'800  K 
and  a  density  of  about  1.42  g/cm^.  Thus  less  energy  is  reieasc-u 
and  a.  lower  initial  pressure  is  created  in  a  given  we!,'ht  ui 
explosive.  In  essence;  it  takes  1  l/5  kg  of  PelletcJ  '  ■  proviae 
an  energy  release  comparable  to  that  from  1  I’.g  of  cn:.,  M.’r. 
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In  addition,  the  detonation  product  gases  are  different  and 
behave  differently  on  expansion.  At  densities  close  to  normal  air 
densities  the  effective  adiabatic  exponent  (7)  is  1.27  for  cast  TNT, 
but  it  is  1.20  for  the  Pelletol.  Since  the  average  pressure  from 
this  gas  in  a  cavity  will  be  proportional  to  the  energy  density 
multiplied  by  the  factor  (7-I),  the  difference  for  the  two  explosives 
is  considerable.  This  effect  is  in  addition  to  the  effect  of  the 
lower  initial  energy  release  of  Pelletol. 

To  insure  a  proper  energy  deposition  on  detonation  and  a  correct 
detonation  velocity,  a  bum  routine  was  provided  in  the  calculations 
which  added  energy  at  an  appropriate  rate  in  the  few  zones  associated 
with  the  detonation  front.  The  position  of  the  front  as  a  function 
of  time  was  fixed  by  the  detonation  velocity,  and  the  rate  of  energy 
deposition  was  dictated  by  the  desired  bum  yield.  The  accompanying 
hydrodynamic  conditions  and  equation  of  state  then  fixed  the  deto¬ 
nation  pressure  and  compression.  Once  the  detonation  had  burned  out, 
the  problem  was  reduced  to  one  governed  by  hydrodynamics  alone. 

Several  calculations  were  made  to  investigate  the  sensitivity 
of  results  to  the  burn  conditions  and  to  the  equation  of  state. 

The  importance  of  the  air  in  the  chamber  (partially  evacuated)  was 
also  Investigated.  Figure  1  compares  the  wall  pressures  for  four 
calculations,  a' 1  for  20  lb  (9  kg)  of  explosive  in  a  6-ft-radlus 
(1.83  meters)  cavity.  The  highest  reflected  pressure  occurs  for  the 
cast  TNT  (Pg  =  1.5)  reaching  570  Ib/ln  (590  bar).  The  rest  are 
for  calculations  using  Pelletol  (p^  »  l.O).  The  sharpest  and  next 
highest  pressure  peak  was  for  Pelletol,  but  instead  of  being  burned 
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with  a  proper  detonation  wave,  all  the  energy  vas  dumped  in 
uniformly  with  all  velocities  zero  initially  and  with  no  initial 
compressions.  In  addition,  this  second  case  had  a  vacuum  in  the 
chamber  (no  air).  A  third  calculation  treated  a  proper  detonation 
and  included  air  in  the  chamber  at  l/20  of  normal  air  density. 

The  last  represented  calculation  (and  the  lowest  peak  reflected 
pressure)  was  also  for  a  detonated  Pelletol  charge  in  0.05  air, 
but  with  a  further  alteration  in  the  equation  of  state  which  changed 
the  behavior  of  the  specific  heat  along  an  adiabat. 

These  four  calculations  (Figure  l)  £-till  do  not  lead  to  an 
unambiguous  comparison  of  the  effect  of  detonation  versus  the  effect 
of  "no-burn"  on  the  resulting  wall  pressure,  and  to  remedy  this,  a 
comparison  calculation  was  run  both  with  and  without  burn  for  the 
case  of  200  lb  (90  kg)  in  a  15-ft  cavity.  A  peak  wall  pressure  of 
59  bar  resulted  for  the  "no-burn"  case  followed  by  a  very  sharp 
drop  to  around  6.9  bar  while  the  case  including  the  detonation 
detail  rose  to  a  flatter  peak  at  around  l4.8  bar.  These  cases 
are  compared  in  Figure  2.  The  peak  pressures  are  obviously  very 
sensitive  to  the  detonation  details,  which  is  in  some  ways  un¬ 
fortunate,  since  more  information  about  the  explosives  is  then 
required  than  is  normally  available. 

The  short-dashed  curve  in  Figure  2  represents  the  time- 
average  of  the  wall  pressure  after  shod  ■i.rrival  and  up  to  each 
instant,  i.e. 


<P) 


(1) 
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Equation  of  state  for  explosives.  The  equations  of  state  for 
both  forms  of  TNT  were  approximate  fits  to  the  adiabatic  data  of 
Tables  3  and  5  of  Jones  and  Miller  [l9^]  (modified  so  as  to  avoid 
negative  slopes  in  the  internal  energy  versus  temperature).  The 
resuj,.tlng  fits  represent  equations  of  state  which  are  incorrect  in 
detail  but  which  produce  the  correct  Initial  and  "final"  states, 
provided  the  influence  of  multiple  shocks  does  not  put  the  gases 
on  adiabats  of  markedly  different  character.  Since  the  larger 
charges  do  not  expand  to  an;;^''thing  like  normal  air  density  in  these 
small  cavities,  and  since  reflected  shocks  may  be  fairly  strong, 
the  proviso  just  injected  is  not  always  observed  in  these  calculations 
to  that  extent,  then,  the  results  at  late  times  may  suffer  from  lack 
of  a  more  precise  equation  of  state  treatment. 

The  following  approximate  fit  to  the  equation  of  state  for 
Pelletol  (P)  was  used,  based  on  the  adiabatic  data  of  Jones  and 
Miller  [l9^,  Table  5]’ 

P  =  ^  (2.89  +  -^)X10"^  (10^°  ergs/cm^)  (2) 

^  .5+V 

where  T  is  in  °K  and  V  is  in  cm^/g. 

E  =  - 'j—  '.v  +  i.655T\lO"^^  (10^°  ergs/g)  (3) 

^  65(>t-T  XIO  ' 

The  later  ammended  form  appears  as 

E*  =  - ^  1.655t\i0"^^  (10^^  ergs/g)  (4) 

P  650V'+T^X10”' 
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where  V  is  talcen  as  unity  for  o  greater  than  1  g/cm^  and  as  equal 
to  the  specific  volume  (V  =  l/p)  for  p  less  than  1.  The  effect  of 
this  change  was  to  reduce  the  effective  gamma  of  the  detonation 
products  at  intermediate  temperatures  and  lov/  densities  (i.e.  in 
the  region  of  1000°  to  j5000°K  and  at  densities  less  than  the  initial 
loading  density).  This  form  of  the  equation  of  state  leaves  the 
detonation  front  and  the  cold  gas  conditions  unchanged  for  the 
detonation  products.  The  effective  (7-I)  or  PV/e  functions 
resulting  from  this  (equation  and  the  earlier  equation  of  state 
fit  (equations  2  and  3)  are  illustrated  in  Figure  us  functions 
of  the  temperature  for  fixed  densities. 

For  cast  TNT  (p^  =  the  fit  used  was  "based  on  Table  3  of 

the  same  reference. 


^T  ~ 


7-73  T 

8500+t\io“^^ 


+  1.2t\io'^'^ 


(10^*^  ergs/g)  (5) 


2.556  $ 


^  (1.15  +  10.96  *  -  288  $ 

(6) 


10^  +  .l4l  $ 
U 


2343  -  7800  +  9237  f^)|  (10^°  ergs/cm^) 


where  $  =  TXIO 

Detonation  characteristics.  The  actual  burn  rate  is  determined 
by  the  equation  of  state  and  loading  density  of  the  explosive .  The 
Chapman-Jouget  conditions  which  apply  at  the  detonation  front  are 
quite  different  for  the  two  TNT  loading  densities  considered. 

'Z 

For  p^  =  1.5;  and  p^  =  1.0  g/cm'^  the  following  values  were  required 
at  the  detonation  front . 
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Table  1. 


Detonation  Characteristics  for  TNT  and  Pelletol 


Physical  Parameter 

Units 

TNT 

Pelletol 

Loading  Density, 

gm/cm^ 

1-5 

1.0 

Detonation  Velocity,  U 

C  J 

m/msec 

6.507 

5.00 

Internal  Energy,  E 

C  J 

10^*^  ergs/g 

5.447 

4.074 

Pressure,  P 

CJ 

10  /  3 

10  ergs/cm 

15.88 

8.224 

Density,  p^j 

g/cm^ 

2.00 

1.421  ^ 

Temperature,  T 

C  J 

"k 

3400 

3800 
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Cavlty  ^ilr.  V/hen  the  cavity  air  Wcas  included  (at  various 
degrees  of  evacuation),  its  equation  of  state  was  determined  in  the 
calculations  by  the  delailed  fit  given  below.  The  equation  of  state 
of  air  is  of  some  importance  here,  since  the  air  around  the  explosive 
is  compressed  and  heated  well  beyond  a  point  where  its  behavior  can 
be  taken  as  that  of  an  ideal  diatomic  gas. 

For  hydrodynamic  calculations  it  is  sufficient  to  consider 
only  the  caloric  relation  between  the  internal  energy  (E),  the 
pressure  (p),  and  the  density  (p)  of  the  gas  in  question.  Such  an 
equation  of  state  for  air  can  be  expressed  in  many  ways;  for  con¬ 
venience  and  simplicity  of  an  analytic  fitting  the  following  form 

has  been  developed  to  cover  all  temperatures  and  a  range  of  densities 

-6 

from  10  times  the  normal  sea-level  density  of  air  to  10  normal 
density.  The  ideal  gas  relation,  E  =  P/p(7-l),  can  be  generalized 
to  any  equilibrium  gas  by  appropriate  definition  of  the  proportionality 
constant  (7-I)  The  following  fit  agrees  with  RAND  and  National 
Bureau  of  Standards  for  air  [Hilsenrath  and  Beckett,  1955;  Gilmore, 
1955]  to  within  5  per  cent  almost  everywhere.  If  we  define  a  quantity 
M-  =  (7+l)/(7-l),  the  pressure  becomes 


P  -  ^ 

■  P-1 

27y+5 

p  =  1  +  +  B  +  C  +  D+E  +  F  +  G+  H+  I 

5y+i 

p  1-0555 

in  which  y  h  (^)(L)  and  ^  =  ^n(p/p^) 


and  the  literal  factors  (A  through  l)  have  the  form 

(1  -  y)u 

(x^  +  x^Ovy"  +  1 


(7) 

(8) 
(9) 


(10) 
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Table  2  gives  the  values  of  the  various  coefficients  for  each  term 
in  p,.  (in  the  table^  the  exponents  of  numerical  coefficients  indicate 
the  povers  of  10,  e.g.^  2.236  means  2.236  x  10"^. )  The  general 
nature  of  this  fit  to  the  equation  of  state  for  air  is  illustrated 
by  Figure  4. 

Computational  method  for  explosion  hydrodynamics.  The  gas 
dynamics  of  these  spherical  explosions  finds  mathematical  description 
in  a  set  of  partial  nonlinear  differential  equations.  The  technique 
used  in  the  numerical  integration  of  these  equations  is  similar  to 
previously  reported  Lagrangean  formulations  in  which  an  artificial 
viscosity  is  used  to  treat  shock  waves  [Von  Neumann  and  Rlchtmyer, 

I95O;  Erode,  1955;  1959] •  An  outline  of  the  numerical  method 
employed  for  the  present  cavity  explosions  is  Included  as  Appendix  A, 
since  much  of  the  detail  is  not  generally  available  elsewhere,  and 
since  an  appreciation  of  the  relative  veracity  of  the  solutions  is 
most  dependent  on  an  understanding  of  the  physical  approximation 


embodied  in  the  mathematical  model. 
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CAVITY  HYDRODYNAMICS 


Results .  The  wall-pressure  histories  were  extracted  from  these 
calculations  to  provide  an  input  to  the  calculations  of  elastic  wave 
propagation.  These  histories  are  illustrated  in  Figures  1^  2,  and 
5  through  9 • 

To  indicate  something  more  of  the  behavior  of  the  blast  in  such 
a  cavity,  some  pressure  density  and  temperature  profiles  for  the  case 
of  20  lb  (9  kg)  in  a  12-ft -diameter  (5.66  meters)  cavity  are 
illustrated  in  Figure?;  10  through  12.  Figure  10  displays  pressure 
as  a  function  of  radius  at  various  times  after  detonation.  In  it  one 
can  follow  the  expansion  from  a  time  near  the  end  of  the  Pelletol 
burning  (t  ~  0.025  msec)  (where  product  gases  have  pressures  near 
100  kb)  through  the  time  when  the  wall  pres8ui*e  reaches  a  maximum 
and  the  reflected  shock  has  begun  its  first  implosion  on  the  origin 
(t  0.78  msec).  At  first  (t  ~  0.0^2  msec)  the  expanding  high- 
explosive  product  gases  ’’snowplow"  the  rarefied  (P^  ~  O.O509  bar)  air. 

But  as  the  expansion  proceeds  (t  ~  0.15^  msec),  the  compressed  air 
begins  to  build  up  an  excess  pressure  ahead  of  the  expanding  explosion 
products.  As  the  restraint  of  the  cavity  wall  is  felt,  the  decelerating 
air  and  explosion  product  gases  pile  up  a  reflecting  shock  (t  ~  O.527  msec) 
and  reach  a  peak  reflected  pressure  at  t  O.7805  msec  of  15  bars 
(188  Ib/in^). 

As  the  pressure  contour  history  continues,  the  first  reflected 
shock  is  seen  to  implode,  reflecting  at  the  origin  to  cause  a  second 
shock  at  the  wall  (at  ~  5  msec),  followed  by  several  rather  weak  shock 


traversals . 
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The  gas  densities  follow  a  similar  history  through  the  detonation 
and  first  reflection  (Figure  ll)  but  lead  eventually  to  a  hotter, 
lower  density  shell  of  air  at  the  wall  most  noticeable  at  late  times. 

The  temperatures  are  initially  somewhat  less  than  4000°K  in  the 
detonation  products  (Figure  12),  cooling  on  expansion,  but  rising 
again  on  being  compressed  by  the  reflected  shock.  Note  that  the 
expanded  gases  drop  to  temperatures  less  than  room  temperature  before 
jumping  to  around  5000°K  through  the  reflected  shock.  The  implosion 
at  the  origin  brings  temperatures  there  of  the  order  of  8000°K. 

This  preferential  heating  leaves  a  temperature  distribution  at  later 
times,  but  before  mixing  and  gravitational  effects  become  important. 

The  air,  compressed  against  the  wall,  is  then  hot  (<^  4500^)  while 
the  detonation  product  gases  show  a  sharp  rise  from  less  than  5000^^ 
to  a  maximum  at  the  origin  close  to  8500°K- 

Since  only  the  wall  pressure  histories  ire  needed  to  define 
the  elastic  wave  propagations,  such  raore  complete  hydrodyniimic 
descriptions  for  the  other  cases  for  which  calculations  were  maue 
have  not  been  included.  Table  3  represents  a  sumnu  vy  of  the  various 
conditions  governing  each  of  the  calculations,  as  well  as  such 
information  as  the  peak  wall  pressures,  time  intervals  between  first 
and  second  shock,  approximate  final  wall  pressure,  length  of  time 
covered  by  the  calculations,  and  appropriate  figure  numbers 
illustrating  the  wall  pressure  histories. 


Pelletol,  P  =  Pelletol^  revised  equation  of  state 
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Comparlson  vith  measured  cavity  vail  pressures.  A  comparison  of 
these  calculated  wal^  pressures  with  the  measured  (gage)  pressures  may 
help  to  put  in  proper  perspective  the  importance  of  the  variables 
mentioned  previously  in  connection  with  choices  of  boundary  and  detona¬ 
tion  conditions^  equations  of  state,  and  inclusion  or  exclusion  of 
rarefied  air.  Figures  15;  l4;  and  15  offer  comparisons  for  the 
9O-,  25O-;  and  450  kg  shots  in  the  9  meter  cavity.  The  data  were 
obtained  only  crudely  from  photographic  reproductions  of  the  instrument 
recordings  and  are  labeled  according  to  the  field  station  designations. 

No  attempt  is  made  here  to  interpret  from  these  the  gage  type  or  location 
(other  than  at  the  cavity  boundary)  or  to  make  any  implications  or 
interpretations  on  the  basis  of  instmimentation  details,  (it  should 
be  reiterated  that  the  data  curves  shown  are  inexact;  for  more  precise 
data  reproductions;  reference  to  the  original  or  direct  copy  tracings 
would  be  advisable.) 

In  Figure  15;  two  data  records  are  indicated  in  comparison  with 
two  calculations.  Of  the  latter,  one  is  for  Pelletol  properly  detonated, 
and  the  other  is  for  Pelletol  without  burn  (as  in  Figure  2).  Beyond 
general  similarities,  one  can  discern  that  the  second  reflection 
occurs  consistently  earlier  (by  about  1  msec)  for  the  calculated  cases. 
The  semilog  scale  of  pressure  should  be  recognized  as  overestimating 
the  differences  at  low  pressures  (where  the  records  are  hard  to  resolve) 
and  as  underestimating  the  differences  at  high  pressures  (where  the  data 
records  seem  to  suffer  from  extraneous  zero-time  transient  signals). 
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Figure  l4  compares  the  calculation  for  2^0  kg  in  the  9  meter 
cavity  with  two  pressure  records,  shoving  again  a  tendency  in  the 
calculation  toward  early  second  pulse  arrival,  or,  more  exactly,  a 
shorter  period  of  oscillation  throughout.  Early  gage  pressures  seem 
to  he  higher  except  in  the  calculated  first  spike.  This  spike  Is  not 
evident  in  the  gage  records.  Aside  from  such  obvious  differences, 
the  calculation  appears  to  provide  the  general  transient  characteristics, 
as  well  as  an  appropriate  average  or  late-time  pressure  value. 

In  Figure  15  comparison  is  made  for  450  kg  in  the  9  meter  cavity. 
Although  the  early-time  comparison  again  shows  higher  measured  pressures 
and  no  spike,  the  period  for  the  calculation  in  this  case  appears  to 
be  longer  than  that  evidenced  by  the  data. 

The  period  discrepancies  may  stem  from  a  combination  of  inexacti¬ 
tudes:  in  both  the  calculations  and  the  actual  shots  the  charge  weights 
and  detonation  conditions  were  not  precise;  for  example,  in  the 
experiments  the  charges  could  not  be  exactly  center-detonated  nor 
could  the  explosive  packing  density  be  exactly  controlled.  Since 
the  packing  density  can  evidently  affect  the  detonation  velocity 
and  energy  release,  the  effective  charge  weight  should  be  considered 
an  approximate  guide  to  the  net  explosion  energy.  Similarly,  in  the 
calculations,  the  detonation  conditions  were  only  approximately 
controlled,  and  detonated  velocities  and  energy  releases  were  neither 
precise  nor  constant.  In  view  of  the  variations  and  vagaries  in  the 
calculations  alone,  the  discrepancies  in  period  and  in  peak  values  of 
pressure  (of  the  order  of  10  to  15  per  cent)  are  not  surprising. 
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But  the  calculated  peak  values  of  the  second  pulse  seem  to  be 
too  low^  even  on  this  loose  scale  of  comparison.  In  the  2^0  kg  case, 
the  calculated  second  peak  is  6.2  bar  while  the  data  curves  range 
from  12.8  to  20.7  bar,  and  in  the  kg  case  the  data  curves  both 
reach  27-6  bar,  while  the  calculation  rises  to  less  than  ll  bar. 
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ELASTIC  RADIATION  CALCULATIONS 


Calculations  of  spherically  symmetric  dilation  waves  resulting 
from  a  time -dependent  pressure  on  the  wall  of  a  spherical  cavity  have 
been  carried  out  with  considerable  generality  by  many  writers 
[Eringen,  1957>  Hopkins,  i960].  The  present  calculations  were 
accomplished  in  a  straightforward  manner  to  describe  the  spherical 
elastic  radiation  field  generated  by  the  theoretical  wall-pressure 
histories  characteristic  of  the  Cowboy  cavity  explosions.  Before 
discussing  these  results  we  shall  indicate  the  method  by  which  they 
were  calculated.  The  basic  equations  are  presented  below  without 
derivation.  A  derivation  of  these  particular  equations  is  given  by 
I^arkin  [1962]. 

Let  a  be  the  radius  of  the  cavity,  P(t)  the  pressure  acting  on 
the  cavity  wall,  and  aF(t)  the  radial  elastic  displacement  at  the 
cavity  wall.  The  function  F(t)  is  given  by 


8ina)(t-T) 


where 


(4) 


^  -  (1-p) 


2 


and 


P  = 


2\x-\ 

2ii+\ 


The  quantities  X.  and  p,  denote  the  Lame  constants  of  the  elastic 
material  and  c  is  the  dllatatlonal  wave  speed. 


(12) 


dT 


The  density  of  the  medium  is  denoted  by  p.  The  radial  strain 


e^(r,t),  therv  is  given  by 


e  F(§)d5 
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The  transverse  strain 


is  given  by 


ac 
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t  -  (r-a)/c 

exp(c§/a)F(?)d§ 


(14) 


From  these  results  one  can  compute  the  radial  elastic  displacement 
u(r,t)  =  rSg,  the  dilatation  +  2eg  =  A(r,t),  and  the  two  stress 
components 


^0  "  ^^r  ^ 

and  ( 15 ) 

a  =  (A.+2(j,)e  + 

r  r  0 

The  particle  velocities  and  accelerations  are  obtained  by 
differentiating  the  fomula  for  u(r,t)  with  respect  to  time. 

All  the  preceding  formulas  were  programmed  for  RAND's 
Johnniac  digital  computer.  Numerical  values  for  the  stresses, 
strains,  dilatation,  particle  displacement,  particle  velocity, 
and  particle  acceleration  were  calculated  as  functions  of  time, 
t,  at  any  radius,  r,  from  the  center  of  the  cavity. 
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In  conjunction  with  the  spherical  wave  code  outlined  above,  a  short 
Fortran  code  has  been  written  for  the  spectral  analysis  of  the 
particle  velocities.  This  program  calculates  the  quantity 


F(u) 


K 

/,N  -iwt,. 
u  ( t ) e  dt 

JOO 

(l6) 


In  addition  to  the  machine  calculations  noted  above,  hand 
calculations  were  made  to  obtain  the  total  energy  radiated  into  the 
elastic  medium. 

Imagine  a  volume  V  of  the  elastic  medium  which  moves  solely  by 
the  deformation  of  the  material.  The  energy  flux  W  radiated  through 
the  surface  S  of  this  volume  is  equal  to  the  rate  at  which  the 
stresses  on  S  do  work  upon  the  material  outside  V  [Thomas,  1961]. 

If  the  time  interval  T  Includes  the  total  pulse  du;ration,  we  can  write 
the  approximate  result  for  spherically  symmetric  motions  so  that 

2 

W  =  4«r  u  a  dt  (17) 

J  U  I* 

o 

where  it  is  assumed  that  r(r^,t)  can  be  replaced  by  r^.  This 
Lagrangean  formulation  does  not  include  any  transport  terms. 

The  elastic  constants  used  in  the  calculations  indicated  above 

are 

Young's  modulus,  E  =  353*^  kb 

Poisson's  ratio,  =  0.25 

Dilatation  wave  speed,  c  =  4.39  x  10^  m/sec 
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After  the  present  calculations  were  started^  Nlcholls  et  al.  [1960] 
published  slightly  different  definitive  values.  The  present  V€LlueS; 
however,  appear  to  be  close  enough  to  those  cited  so  that  no  further 
calculations  are  necessary. 

Calculations  were  carried  out  for  cavity  radii  of  I.83  and 
4.57  meters  corresponding  to  ten  of  the  twelve  hydrodynamic  calculations. 
The  two  cases  for  explosions  of  cast  TNT  were  excluded  from  further 
study,  as  they  were  calculated  for  time  intervals  too  short  to  be  of 
present  use.  For  each  theoretical  P(t)  from  these  Pelletol  explosions, 
elastic  wave  calculations  were  run  for  three  radii  which  were  chosen 
to  approximate  closely  the  experimental  gage  locations.  Thus,  for 
the  1.83  meter  cavity,  r  values  were  chosen  to  be  at  3-05;  6.3I,  and 
11  meters  from  the  center  of  the  cavity.  For  the  4.57  meter  cavity, 
chosen  r  values  were  J-62,  15*4,  and  24.4  meters. 
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THEQRETICAL  ./VND  EXPERIMENTAL  RADIATION  DATA 

As  explained  above,  the  hydrodynamic  conditions  within  the 
cavity  are  markedly  Influenced  by  the  nature  of  the  burning  process 
and  by  the  equation  of  state  for  the  explosive  products.  These 
effects  are  directly  mirrored  by  the  response  of  the  elastic  medium, 
as  evidenced  by  the  large  variations  in  radiated  energy  of  Table  ^ 
for  the  same  charge  weight  but  for  the  different  equations  of  state 
for  the  explosive.  Of  course,  the  calculated  values  of  the  radiated 
energy  depend  upon  the  time  Interval  over  which  the  calculation  is 
extended.  Even  for  those  cases  of  longest  duration  it  is  not  strictly 
true  that  the  entire  elastic  pulse  could  be  used  in  (17)*  On  the 
other  hand,  it  was  found  that  the  calculations  of  shortest  duration 
contain  about  90  per  cent  of  the  radiated  energy.  The  calculations 
of  greatest  duration  appear  to  give  energy  values  which  are  accurate 
to  within  1  per  cent. 

Figures  l6  and  17  compare  particle  velocities  in  the  salt  after 
detonation  of  90-  and  2^0-  kg  Pelletol  charges.  As  noted  previously, 
these  theoretical  curves  result  from  using  the  values  of  wall-pressure 
versus  time  from  Figures  2  and  8  as  input  data  for  the  elastic  wave 
program.  Figures  l8  and  19  compare  particle  velocities  and  particle 
displacements  at  two  distances  from  the  center  of  the  4.57  meter 
cavity  corresponding  to  the  theoretical  P(t)  curve  of  Figure  9  for  a 
charge  weight  of  450  kg.  The  two  radial  distances  are  7*62  and  24.4 


meters,  respectively. 


The  over-all  agreement  between  experiment  and  theory  shown  by 
these  four  illustrations  is  thought  to  be  quite  satisfactory.  The 
differences  between  theoretical  and  experimental  values  of  peak  wall 
pressure  from  the  cavity  hydrodynamics  appear  here  as  corresponding 
differences  in  peak  particle  velocities.  The  somewhat  better  agree¬ 
ment  between  experiment  and  theory  for  the  magnitudes  of  subsequent 
peaks  as  well  as  the  improved  correlations  of  the  time  intervals 
between  peaks  is  also  noted  in  the  wall-pressure  curve.  In  Figure 
19  the  difference  between  the  experimental  and  theoretical  arrival 
times  from  the  instant  of  detonation  until  reception  of  the  first 
signal  is  about  1  msec  at  r  =  2k. h  meters.  Figure  I8  and  the  wall- 
pressure  curves  of  Figure  9  indicate  that  this  difference  should  be 
about  1/3  msec.  Inspection  of  the  "digitized"  velocity  data 
[Adams  and  Allen,  I961]  at  r  =  24.4  for  shot  10  of  the  Cowboy  tests 
tests  indicates  agreement  with  the  1/3  msec  figure. 

The  experimental  data  of  Figures  16  through  I9  of  this  paper 
have  been  taken  from  copies  of  the  oscillograph  recordings  which 
were  made  from  the  original  taped  records.  Evidently  the  detonation 
time  marked  on  the  particular  record  at  r  =  24.4  meters  has  been 
inadvertently  misplaced.  Had  we  made  use  of  this  information  in 
Figui*e  19,  the  agreement  between  experiment  and  theory  would  have 
been  greatly  improved.  However,  we  have  applied  no  corrections  to 
the  experimental  data.  The  experimental  displacement  data  of 
Figures  I8  and  19  have  been  read  from  records  that  were  obtained 
from  the  velocity  tapes  by  a  process  of  electronic  integration. 
[Murphey,  19^1 ] •  Comparison  of  experimental  and  theoretical  dis¬ 
placement  curves  indicates  the  possibility  of  some  drift  in  the 


experimental  records.  It  is  to  be  noted  that  Herbst  et  al.[l96l] 


present  an  experimental  displacement  curve  in  Figure  19  of  their 
paper  which  corresponds  to  the  conditions  of  our  Figure  19.  Their 
experimental  curve  results  from  the  Integration  of  the  velocity  record 
(presumably  by  direct  numerical  con^jutation)  and  does  not  show  any 
drift.  Use  of  their  curve  in  Figure  I9,  instead  of  that  read  from 
the  electronically  integrated  record,  would  provide  excellent  agree¬ 
ment  between  experiment  and  theory  for  particle  displacements  at 
r  =  24.4  meters. 

Figure  20  compares  velocity  spectrums  calculated  from  the 
experimental  and  theoretical  velocity  data  of  Figures  I8  and  I9. 

At  high  frequencies  these  spectrums  show  fair  agreement.  The 
pulsation  frequency  of  the  gas  in  the  cavity  is  clearly  shown  by  the 
peak  in  these  curves  between  100  and  200  cps.  The  spectrums 
calculated  from  the  measured  velocities  have  considerably  greater  low- 
frequency  content  than  the  theoretical  spectrums.  Investigation  of 
this  difference  indicated  that  it  was  partly  due  to  the  fact  that 
the  experimental  velocity  curves  of  Figures  18  and  19  have  been 
drawn  with  rather  broad  initial  peaks.  Actually,  as  already  dis¬ 
cussed  by  Murphey  [I961],  the  velocity  records  of  the  decoupled 
shots  showed  a  high-frequency  component  due  to  ringing  of  the 
cannisters  containing  the  gages.  Although  this  high  frequency  on 
the  recordings  is  of  fairly  small  amplitude,  it  obscures  the  width 
of  the  initial  velocity  pulse.  In  the  process  of  trying  to  smooth 
out  these  oscillations  in  drawing  the  experimental  curves,  the 


width  of  the  initial  pulse  tends  to  be  overestimated. 
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For  example,  when  a  somewhat  narrower  initial  pulse  was  obtained  from 
a  rereading  of  the  data  of  Figure  l8,  the  low-frequency  part  oF  the 
experimental  curve  of  Figure  20  was  lowered  by  about  a  factor  of  2. 

No  attempt  has  been  made  to  improve  the  agreement  between  experiment 
and  theory  by  applying  this  factor  because  of  the  uncertainties 
associated  with  such  procedures. 

Another  cause  for  disagreement  between  the  experimental  and 
theoretical  spectnoms  resides  with  the  truncation  of  the  records 
which  were  analyzed.  For  consistency,  all  curves  of  Figure  20  were 
calculated  for  20  msec  after  the  arrival  time  for  the  velocity 
function  u  in  (l6).  Outside  this  20  msec  interval  the  velocity  is 
assumed  to  vanish.  On  the  other  hand,  when  the  record  is  truncated 
at  some  arbitrary  time  T  later  than  the  arrival  time,  (l6)  shows 
for  low  frequencies  (w  =  0,  to  be  exact)  that  F(u))  is  nearly  equal 
to  the  particle  displacement  u(t).  If  the  interval  of  time  used  in 
(l6)  is  short  enough,  u(t)  can  show  considerable  variation,  depending 
upon  the  choice  of  T.  The  largest  variation  in  F(w)  resulting  from 
this  cause  was  found  to  be  within  50  P^r  cent  at  u)  =  0. 

The  preceding  comparisons  have  been  made  for  experimental 
conditions  in  which  the  assumptions  of  elasticity  are  reasonably  well 
satisfied.  For  the  larger-yield  explosions  in  the  smaller  sphere, 
the  salt  shows  an  anelastic  behavior,  which  can  be  illustrated  clearly 
by  a  comparison  of  experimental  and  theoretical  results.  For  this 
comi)arison  it  seems  sufficient  to  compare  peak  values  of  the  particle 
velocity  and  of  the  particle  displapement  for  each  case.  The  various 
values  are  entered  in  Table  4.  Nearly  all  of  the  experimental  values 
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given  in  Table  4  have  been  reported  by  Murphey  [1961]  in  his  Tables 
3  and  4.~  The  first  thing  indicated  by  Table  4  is  that  the  peak 
velocities  and  displacements  deduced  from  the  Pn  calculations  show 
considerable  differences  from  the  experimental  quantities,  while  the 
values  calculated  from  the  Pb  and  P'b  wall  pressure  functions  are  in 
fair  or  good  agreement  with  experiment.  A  second  thing  revealed  by 
these  data  is  that,  for  those  cases  in  which  elasticity  of  the  salt 
is  a  good  approximation,  the  agreement  between  experiment  and  theory 
seems  to  be  best  at  those  radii  closest  to  the  center  of  the  cavities. 
Comparison  of  trends  at  greater  ranges  appears  to  Indicate  that,  even 
for  small  strains,  the  salt  has  an  element  of  dissipation  that  is  not 
accounted  for  by  the  theory  .  Finally,  the  effects  of  plasticity  or 
other  anelastic  behavior  are  shown  dramatically  by  the  results  from 
the  "overdriven"  cavity  of  I.83  meter  radius  with  the  9OO  kg  charge. 

Comparisons  between  experimental  and  theoretical  data  for  other 
physical  quantities  are  possible.  The  preceding  comparisons,  however, 
seem  to  be  the  most  useful  ones  at  this  time.  For  example,  the  particle 
acceleration  data  from  the  tests  shov'  exceedingly  large  variations  from 
the  theoretical  values,  and  the  body  of  these  measurements  does  not 
appear  to  be  internally  consistent .  Similarly,  the  data  from  the 
"earth-pressure"  gages  (whatever  it  is  they  measure)  show  great  scatter 
among  themselves,  and  these  pressures  show  no  resemblance  to  calculated 
values  of  a^,  or  ~  ^0 ’  Unfortunately,  e^(t)  data  from  the  dynamic 
strain  gages  [Nlcholls  et  al.,  I960]  were  obtained  for  r  values  that 
were  considerably  larger  than  any  radius  used  in  present  calculations. 

^In  Table  4  of  Murphey  [1961]  shot  5  should  be  designated  as 
shot  2  and  shot  4  should  be  designated  as  shot  3* 


Table  4.  Experimental  and  Calculated  Peak  Particle  Velocities  and  Displacement 
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For  this  reason^  no  direct  comparison  between  theoretical  and 
measured  radial  strains  is  possible.  Even  so,  the  rise  times 
for  the  experimental  strains  appear  to  be  considerably  longer 
than  those  that  would  be  found  from  elastic  calculations.  This 
finding  supports  the  evidence  from  Table  4  concerning  the  effect 
of  dissipation  at  small  strains  in  the  salt. 


CORRELATIONS  FOR  CLOSE-IN  EFFECTS 


The  complex  nature  of  the  sequence  of  events  inside  the  cavity 
makes  it  difficult  to  find  useful  ways  of  correlating  characteristic 
features  of  the  waves  in  the  salt.  However,  if  such  a  systematization 
of  the  data  can  he  achieved,  the  theoretical  findings,  which  are  free 
of  extraneous  influences,  should  offer  the  greater  chance  of  success. 
Even  if  such  correlations  are  found,  it  is  important  to  remember  that 
they  will  be  based  upon  the  approximation  that  salt  is  an  elastic  solid- 
In  situations  for  which  the  preceding  results  have  shown  that  this  is 
not  actually  a  good  approximation  for  salt,  the  fiction  of  elasticity 
is  maintained  below  because  of  its  theoretical  interest. 

Table  3  permits  the  comparison  of  the  energy  released  by  the 
explosions  to  that  radiated  into  the  salt.  A  significant  feature 
of  these  results  is  that  the  radiated  energy  W  is  always  a  very  small 
fraction  of  1  per  cent  of  U,  the  explosive  energy  released.  To  obtain 
a  graphic  idea  of  this  result,  one  can  represent  W  as  a  function  of  U 
in  a  log-log  plot  to  find  that  the  radius  of  the  cavity,  a,  and  the 
nature  of  the  burn  calculation  for  the  explosive  affect  the  result. 
Otherwise  the  radiated  energy  appears  to  be  proportional  to  a  constant 
power  of  the  explosive  energy.  It  is  found  that  the  radiated  energy 
decreases  as  the  cavity  radius  increases,  so  that  one  is  Led  to  plot 
W  as  a  function  of  the  ratio  u/a.  The  result  of  this  plot  is  shown 
in  Figure  21  where  now  only  data  corresponding  to  explosions  with  or 
without  detonation  calculations  lie  on  separate  lines.  Actually  the 
data  for  the  modified  equation  of  state  (P'  )  He  on  or  below  the  solid 
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line  in  Figure  21,  while  the  data  for  the  equation  of  state  designated 
by  P  are  above  the  line,  so  that  even  in  this  case  some  differences  in 
the  details  of  the  burn  are  evident.  However,  as  noted  previously, 
the  limited  accuracy  of  the  calculations  for  the  radiated  energy  W 
does  not  permit  a  definite  distinction  between  these  two  cases. 

The  curves  of  Figure  21  suggest  that  for  detonated  Pelletol  we  have, 
approximately, 

Vf  a  0.6  X  10  ^  (u/a)^*^^  ergs  (l8) 

and  for  the  (Pn)  calculations,  which  exclude  detonation,  we  have, 
approximately, 

W  =  2  X  10  (u/a)^’^^  ergs  (19) 

for  which  u/ a  is  expressed  in  dynes. 

Because  the  elastic  constants  of  the  solid  were  not  varied  in 
the  present  theoretical  calculations,  it  is  not  surprising  that  the 
preceding  correlation  shows  no  explicit  dependence  upon  these  constants. 
On  the  other  hand,  when  the  elastic  constants  for  salt  were  used  in  an 
effort  to  obtain  dimensionless  correlations  in  accordance  with  the  well- 
known  laws  of  similitude  [Parkin,  1958] >  no  correlation  of  these  results 
was  found.  For  this  reason,  calculations  for  solids  with  different 
elastic  constants,  but  for  the  same  hydrodynamics  inside  the  cavity, 
might  be  of  considerable  help  in  obtaining  a  more  general  result. 

Although  the  radiated  energy  W  appears  to  be  simply  related  to 
the  energy  released  by  the  explosive  U,  the  essentially  different 
magnitudes  of  these  two  quantities  suggest  that  two  characteristic 
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lengths  might  be  significant  for  the  dynamic  response  of  spherical 
cavities.  Thus,  for  effects  associated  with  the  explosion  Itself, 
we  define  a  characteristic  length, 

a  =  (u/e)^'^^  (20) 

and  for  events  out  in  the  elastic  medium  we  define  another  characteristic 
length, 

P  (w/e)^/^  (21) 

(The  length  p  is  not  to  be  confused  with  the  parameter  p  defined 
previously  in  connection  with  equations  11  and  12.) 

The  Young's  modulus  E  is  taken  as  a  characteristic  pressure  for  both 
a  and  p  even  though  it  might  seem  more  appropriate  in  the  case  of  a 
to  use  a  pressure  characteristic  of  the  explosive.  Again  we  note  that 
the  elastic  constants  of  the  solid  have  not  been  varied  in  these 
calculations,  so  that  the  use  of  E  in  (20)  and  (21)  secur-es  only 
dimensional  consistency. 

The  characteristic  lengths  a  and  p  are  useful  for  graphical 
presentations  of  the  peak  values  of  calculated  quantities.  For 
example.  Figure  22  is  a  dimensionless  plot  of  peak  displacements 
versus  the  distance  from  the  center  of  the  cavity.  The  open  points 
denote  results  calculated  from  the  P'b  equation  of  state,  the  solid 
point  symbols  indicate  results  from  the  Pb  calculations,  and  symbols 
which  are  half  solid  and  half  open  designate  Pn  results.  This 
symbolism  has  been  employed  for  Figures  23,  24,  and  25  also.  Figure  22 
is  the  only  instance  in  which  it  was  found  possible  to  achieve  a 
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reasonably  close  correlation  of  the  calculated  results  in  terms  of 
dimensionless  quantities.  For  example^  in  Figure  it  vas  necessary 
to  plot  the  product  of  peak  particle  Mach  number  u  /c  and  3  ^  against 
the  scaled  distance  r/a.  Similarly,  in  Figure  2k  the  best  over-all 

correlation  is  obtained  when  the  product  of  dimensionleLS  acceleration 

2  -2 

u  3/c  and  3  Is  plotted  versus  scaled  range.  In  this  last  case 

t  O 

the  correlation  is  not  as  good  as  that  obtained  in  the  two  preceding 
examples.  In  Figure  25  the  peak  radial  stresses,  including  the  peak 
cavity  wall  pressures,  are  shown  as  a  function  of  scaled  range. 

The  most  efficacious  combination  of  parameters  for  this  correlation 
was  found  to  be  a ^E3. 

In  view  of  the  difficulties  involved  in  the  scaling  of  single 
peak  values  of  the  various  quantities  studied,  it  is  thought  that 
the  results  given  above  constitute  fairly  good  correlations.  The 
fact  that  the  points  from  the  Pn  calculations  fall  among  the  points 
from  the  Pb  and  P'b  calculations  is  most  encouraging.  Because  the 
system  contains  two  characteristic  lengths,  the  task  of  finding  such 
trends  is  rather  tedious,  and  we  have  not  explored  all  possible 
combinations.  Furthermore,  considerations  of  similitude  [Park in , 

1958]  indicate  that  dimensionless  quantities  such  as  u/3, 
u^^3/c  ,  and  CT are  the  appropriate  parameters  for  such  correlations, 
and  these  quantities  have  been  employed  even  though  the  elastic 
constants  have  not  been  varied  in  the  present  study.  For  the  same 
reason,  these  parameters  have  been  retained  even  when  it  was  not 
possible  to  obtain  correlations  with  respect  to  these  dimensionless 


quantities . 
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From  the  foregoing  we  see  that  if  the  dimensionless  qxiantities 

2  n 

u/p,  u  /c;  and  u  p/c  are  multiplied  by  the  quantity  0  (where  n  = 

0,  1,  or  2  for  the  peak  displacements,  velocities,  or  accelerations, 

respectively),  these  quantities  appear  to  be  inversely  proportional 

to  the  dimensionless  range  r/a.  For  the  dimensionless  radial  stress, 

the  appropriate  value  for  n  is  unity. 

The  only  data  available  to  test  the  utility  of  the  present 
correlation  schemes  consist  of  the  measured  peak  particle  displacements 
and  velocities  of  Table  4.  In  the  absence  of  experimental  measurements 
for  the  energy  radiated  into  the  salt,  we  have  calculated  a  value  for 
each  shot  of  the  Cowboy  tests  as  follows:  First,  the  explosive  energy 
U  for  each  charge  was  calculated  by  using  the  value  of  4.074  x  10^^ 
ergs/g  for  the  energy  density  of  Pellet ol.  Then  (l8)  was  used  to 
calculate  the  radiated  energy  W.  The  appropriate  values  of  a  and  3 
for  each  shot  follow  from  (20)  and  (2l).  With  the  help  of  these  data 
and  the  sound  speed  in  the  salt,  the  experimental  data  of  Table  4  can 
be  characterized  in  terms  of  the  parameters  u/p,  u  /c3^  and  r/a- 
Plots  of  these  quantities  appear  in  Figures  26  and  27.  The  correlation 
of  these  data  appears  to  be  quite  good.  It  suggests  approximate 
power-law  relationships  of  the  form 

u/3  =  2.l(r/a)"^'®  (22) 


and 

u  /c3  =  0.l86(r/a)"^‘^ 
0 


per  meter 


(25) 
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It  is  interesting  to  note  that  the  experimental  data  show  a 
more  pronounced  attenuation  of  peak  scaled  particle  displacement 
and  velocity  with  increasing  scaled  range  than  do  the  corresponding 
theoretical  results.  The  theoretical  values  appear  to  be  approximated 
simply  by  an  inverse  proportionality.  Two  factors  that  are  probably 
responsible  for  this  difference  are  those  of  plasticity  and  of  other 
dissipation  in  the  salt  which  occurs  even  at  small  strains.  The 
first  of  these  factors  is  most  important  at  small  scaled  ranges. 

For  a  given  explosive  energy^  plasticity  in  the  salt  peimiits  larger 
particle  displacements  and  velocities  than  would  be  obtained  if  the 
salt  were  truly  elastic.  The  dissipation  at  small  strains  is  most 
effective  at  large  scaled  ranges  where  it  tends  to  diminish  the 
magnitude  of  the  peak  displacements  and  velocities  more  strongly 
than  it  would  in  an  elastic  material.  Of  course^  the  calculation  of 
the  characteristic  length  6  has  been  carried  out  under  the  supposition 
that  the  salt  is  an  elastic  material.  To  the  extent  that  the  salt 
departs  from  this  idealized  behavior,  the  results  of  Figure  l8  are 
inaccurate.  Moreover,  the  use  of  theoretical  values  of  the  radiated 
energy  instead  of  measured  ones  is  also  a  source  of  error. 
Nevertheless,  Figures  26  and  27  serve  to  illustrate,  in  a  fairly 
quantitative  way,  how  wave  propagation  in  the  hau.ite  of  the  Cowboy 
tests  differs  from  that  which  would  take  place  in  an  elastic  solid. 


CONCLUSIONS 


For  the  explosion  phenomena  within  the  cavities  the  following 
features  were  investigated  and  found  to  have  significant  influence 
on  cavity  wall  pressures;  (l)  loading  density  of  explosive^ 

(2)  equation  of  state  of  explosive  products  and  of  air,  (5)  the 
presence  of  air  In  the  cavity,  and  (4)  detonation  details  of  explosion 
initiation  and  burn. 

To  a  similar  extent,  the  control  of  these  same  features  must 
inevitably  influence  the  actual  experiments.  Although  great  care 
went  into  the  experiments  to  control  and  to  record  as  many  features 
of  the  explosions  as  practical,  the  very  nature  of  explosives  and 
detonations  is  in  conflict  with  notions  of  precision  and  exact 
reproducibility,  to  the  end  that  when  coupled  with  the  necessary 
instrumentation,  uncertainties  and  comparisons  of  wall  pressure 
measurements  with  calculated  wall  pressures  for  any  one  shot  are 
not  easily  interpretable.  In  addition,  several  other  "non-ideal” 
effects  were  of  necessity  present  in  the  experiments:  The  initiations 
were  only  approximately  central,  and  not  entirely  spherically  uniform. 
The  existence  of  supporting  platforms  and  tripods  together  with 
packing  material  in  the  top  of  the  plastic  explosive  containers  may 
have  added  to  the  lack  of  symn^try. 

The  influences  of  a  few  important  features  of  such  an  explosion 
have  been  investigated  here,  and  it  is  hoped  that  an  appreciation  of 
these  factors  may  aid  both  in  the  interpretation  of  these  experiments 
and  in  the  planning  of  futxire  explosives  tests. 
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The  calculations  of  elastic  waves  presented  here  complete  our 
attempt  to  give  a  wholly  theoretical  set  of  calculations  for  the 
cavity  explosion  phase  of  the  Cowboy  program.  In  view  of  the  highly 
idealized  theoretical  model  and  of  the  uncertainties  concerning  the 
properties  of  the  explosives  used  in  the  tests,  the  over-all  agree¬ 
ment  between  close-in  experimental  measurements  and  theoretical 
calculations  is  thought  to  be  excellent.  Of  course,  those  explosions 
which  seriously  overstressed  the  cavities  must  be  excepted.  These 
overdriven  cavity  examples  when  compared  with  the  corresponding 
elastic  calculations,  give  a  good  indication  of  the  close-in  effects 
of  non-elastic  response.  Even  in  those  cases  for  which  the 
assumption  of  elasticity  is  a  good  approximation,  comparison  of 
experimental  and  theoretical  trends  has  indicated  some  dissipative 
effects  in  the  salt  which  are  not  present  in  an  elastic  material. 
Finally,  the  theoretical  calculations  have  indicated  the  existence 
of  some  correlations  between  (l)  explosive  energy  and  radiated 
elastic  energy  and  (2)  closs-ln  distances  and  the  peak  elastic 
displacements,  velocities,  accelerations,  and  radial  stresses. 

These  correlations  could  be  of  considerable  use  for  predictions  of 
these  quantities  for  other  cavity  explosions  in  elastic  materials. 
Before  this  potential  usefulness  can  be  reliably  exploited,  however, 
additional  calculations  should  be  carried  out  for  the  same  and 
other  explosive  sources  in  other  elastic  materials  in  order  to 
generalize  the  present  findings. 

Furthermore,  we  note  that  the  correlations  of  experimental 
dimensionless  peak  particle  displacements  and  scaled  peak  particle 
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Mach  nionibers  with  scaled  range^  although  most  encouraging,  would  be 
much  more  meaningfuJ.  if  they  were  based  upon  experimental  values  of 
the  radiated  energy.  The  fundamental  importance  of  the  radiated 
energy  should  make  its  experimental  determination  an  essential  part 
of  future  cavity  explosion  experiments.  In  such  experiments,  in 
elastic  or  anelastic  solids,  the  fundamental  quantity  is  the  energy 
radiated  from  the  surface  of  the  cavity.  Therefore,  if  one  can 
measure  the  velocity  of  the  cavity  wall  (or  a  particle  velocity  in 
the  solid  sufficiently  near  the  wall)  as  a  function  of  time,  this 
datum  together  with  the  wall  pressure  history  will  permit  a  determina¬ 
tion  of  the  radiated  energy.  We  hope  that  a  feasible  experimental 
technique  for  this  determination  will  be  found. 


Appendix .  COMPUTATIONAL  METHOD 


The  problem  presented  by  explosions  lies  in  the  solution  of  a 
set  of  nonlinear^  partial  differential  equations  representing  the 
conservation  of  masS;  momentum^  and  energy.  These  conservation  laws 
may  be  expressed  mathematically  in  several  ways,  but  they  ai’e 
generally  formulated  in  terms  of  either  Eulerian  or  Lagrangean 
coordinates.  The  Eulerian  form  is  an  expression  of  the  conservation 
lavs  as  viewed  from  coordinate  systems  fixed  in  space^  and  the 
Lagrangean  form  is  an  expression  of  the  same  conservations  in  terms 
of  a  fixed  set  of  masses  or  gas  particles,  so  that  a  solution  in  the 
Eulerian  case  represents  the  history  of  the  blast  wave  at  a  fixed 
point  and  a  solution  in  the  Lagrangean  system  describes  the  experience 
of  a  particle  (or  a  volume  of  gas)  as  it  moves  about-  The  latter 
description  is  particularly  fitting  to  this  problem  since  the  burning 
of  an  explosive  is  most  easily  described  in  terms  of  the  histoiy  of 
each  burning  mass. 

The  set  of  equations  below  is  one  representation  of  these 

conservation  lavs  in  a  Lagrangean  system  for  a  spherically  symmetric 

explosion.  In  these  equations  the  Lagrangean  variable  or  mass 

coordinate  is  characterized  by  the  mass  per  steradian;  i.e., 

m  =  p  R'^/3  where  p,  is  the  initial  density  of  the  gas  and  R  is  its 
1  o  1  0 

initial  position. 

i  =  V 

p  3  dm 


conservation  of  mass;  (Al) 


du 

dt 


Is  (P  -  Q)  , 


conservation  of  momentum;  (A2) 


-40- 


9t 

E  =  E(T 
or 

P  =  P(E 


particle  velocity.  (A5) 

In  these  equations  the  symbols  are  defined  as  follows: 

V,  specific  volume 
R;  radial  distance 
u,  particle  velocity 
P,  pressure 

artificial  viscosity  pressure 
E,  specific  internal  energy 
T,  temperature 

Dj  specific  energy  source  (or  loss)  rate 
m,  mass/ steradian 
t,  time. 

In  solving  these  equations  for  a  specific  set  of  initial  and 
boundary  conditions  by  numerical  means,  the  usual  approach  is  to 
arrange  a  set  of  difference  equations  which  approximate  these 
differential  equations  over  a  set  of  finite  mass  elements.  In  such 
a  system,  then,  derivatives  become  differences,  and  a  solution  is 
accomplished  by  the  repetition  of  a  prescribed  series  of  numerical 
steps  which  solve  the  difference  equations  for  new  values  of  the 
parameters  w^hlch  would  occur  after  a  small  time  increment. 


(p  +  Q)  f  *  D  , 


conservation  of  energy;  (A3) 


,V),  P  =  P(T,V) 


equations  of  state;  (a4) 


,V),  T  =  T(E,V) 
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To  formalize  the  differencing  procedure,  consider  the  t:rans- 
formation  of  the  above  differential  equations  to  another  set  in  which 
the  independent  variables  are  a  dimensionless  Lagrangean  coordinate 
J  and  a  dimensionless  time  n,  and  in  which  the  masses  are  represented 
by  Integer  values  of  j  and  the  time  by  integer  values  of  n. 

The  Jacobian  for  the  Lagrangean  coordinate  system  then  follows 
from  the  differential  equation  representing  conservation  of  mass: 

“  3V(J,0)  ^  “  3V(J,n)  ^  (a6) 


The  othex"  conservation  equations  become 


.  .  r2  ^  (P  ^  Q) 


(momentum)  (A7) 


dn 


(energy)  (a8) 


and 


dR  /dt. 


-1 


(A9) 


In  these  equations  the  dependence  of  the  variables  on  j  and  n  has 
been  suppressed  for  simplicity. 

The  artificial  viscosity  may  be  variously  defined,  but  the 
following  is  a  usual  form: 


Q  =  -  (-) 
^  V  ''dj^ 


2 


if  <  0 


du 


i  0 


(AlO) 


Q  =  0 


if 
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In  which  the  value  of  C  determines  the  amount  of  spreading  of  shocks. 
When  C  is  unity,  the  spread  is  about  three  zones,  and  when  C  is  six 
the  spread  is  around  six  zone«,  the  constant  being  proportional  to 
the  square  of  the  spread  [Von  Neumann,  1950;  Erode,  1955] • 

In  some  unusual  circumstances  this  form  fails  to  provide  the 
necessary  shock  spreading.  Such  a  case  can  arise  in  some  problems 
involving  imploding  shock  waves. 

A  more  rigorous  form  for  the  artificial  viscosity,  which  is 
generally  equivalent  to  the  form  above  (equation  AlO),  is  the  following: 


Q 


dV 

0  for  ^  ^ 
dn 


The  formation  of  a  set  of  centered  difference  equations  is  now 
a  matter  of  assigning  integer  values  to  the  independent  variables 
J  and  n  and  either  defining  the  dependent  variables  at  these  points 
or  at  midpoints  (at  half-integer  points)  in  J  or  n.  It  is  conventional 
to  express  the  J  and  n  values  associated  with  a  given  parameter  by 
subscripts  and  superscripts;  thus  the  radius  at  a  mass  point  J  and  at 
a  time  t  after  n  time  increments  At(n)  =  At^  would  appear  as  R^. 

n-f^ 

The  equations  for  the  nth  time  step  At  ®  over  a  set  of  mass 
elements  follow. 


n+i  n-^ 


Amj 


2  ^ 


n  pH  n-|  n-^ 
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for  which 
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Q"i  --  0  for  u""^  >  u""^ 
J-^  t)  J 


E"")  =  E"  -  (i  p^:i  4  p-.i "  ^ 


In  these  equations  it  Is  worth  noting  that  the  quantities  are  either 
defined  or  averaged  In  such  a  way  that  first-order  errors  in 
Am  are  avoided;  l.e.,  the  difference  equations  are  centered.  The 
actual  progression  of  the  solution  can  better  he  illustrated  by 
reference  to  a  space-time  grid,  as  below,  showing  the  points  at 
which  the  various  parameters  are  defined. 


j-1  J'-s  J 
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For  illustration,  assume  that  the  nth  time  step  has  been  completed, 
that,  therefore,  all  the  values  at  times  less  than  n  are  known, 
and  that  the  next  step  is  to  determine  the  values  at  n+|  and  n+1 
times.  The  difference  equations  have  already  been  written  in  the 
order  in  which  they  are  computed.  The  first  equation  defines  a  new 
velocity  at  time  n+^,  at  point  j,  in  terns  of  the  old  velocity  at 
that  j  point  at  time  n-^,  and  the  radius  at  that  j  point  at  time  n, 
together  with  the  pressure  gradients  at  that  point  and  at  time  n  as 
represented  by  the  differences  in  the  pressures  at  J  points  one  half 
a  unit  on  either  side  of  J.  From  this  new  velocity  a  new  radius  is 
computed  in  a  straightfoivard  manner,  using  the  definition  of  the 
particle  velocity.  These  new  radii,  in  turn,  are  used  to  compute 
a  new  specific  volume,  and  that  new  volume  together  with  the  new 
values  of  the  velocity  can  be  used  to  compute  a  new  artificial 
viscosity.  The  artificial  viscosity  is  nearly  of  second  order  in 
the  time,  and  its  lack  of  symmetry  in  the  velocity  equation  does  not 
disturb  the  degree  of  accuracy  of  these  otherwise  centered  difference 
equations.  The  new  volume,  together  with  the  new  artificial  viscosity, 
can  then  be  used  in  the  energy  equation  to  determine  a  new  energy  and 
a  new  pressure.  If  the  equation  of  state  is  not  simple,  it  is  then 
necessary  to  arrive  at  the  new  energy  and  pressure  by  some  iterative 
method.  The  equation  of  state  represents  a  second  relation  between 
these  remaining  two  unknowns,  i.e.,  the  new  energy  and  the  new 
pressure.  With  a  satisfactory  solution  for  all  of  these  parameters 
at  a  point  j  one  can  then  advance  to  the  next  j  point  and  repeat 
the  same  procedure. 
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Of  course,  if  the  prohlem  involves  other  mechanisms  of  energy 
transfer,  or  sources  or  sinks  of  energy,  the  term  D  in  the  energy 
equation  becomes  useful.  For  instance,  the  D  term  can  be  used 
satisfactorily  in  detonation  wave  calculations.  A  number  of  forms 
to  represent  the  buniing  of  high  explosives  or  gases  in  a  detonation 
front  are  possible;  one  fairly  crude  but  workable  form  is  the 
following: 


I,n.l  "c/* 
-  t 


n+i 


Cj 


(A18) 


S  AR 


Cj 


in  which  E  .  is  the  Chapman-Jouget  energy,  l.e.,  the  energy  generated 
C  j 

in  the  detonation,  U  .  is  the  detonation  velocity,  and  S  is  the  spread 

0  J 

or  thickness  of  the  detonation  front  in  keeping  with  the  shock  spreading 
characteristic  of  the  artificial  viscosity  technique. 

A  great  many  cautions  resulting  from  experience  and  from  the 
details  of  the  handling  of  the  numerical  procedure  might  well  be 
included  in  a  complete  description  of  a  numerical  method,  since  many 
special  considerations  are  necessary  before  a  numerical  program  becomes 
truly  workable.  Some  of  these  factors  are  concerned  with  the  stability 
of  the  numerical  procedure,  attention  to  boundary  conditions  and  initial 
conditions,  and  internal  checks  on  the  accuracy  of  the  computation. 

Such  detail,  however,  is  not  considered  appropriate  to  the  completeness 


of  this  paper. 


Fig.  1  -  Comparison  of  cavity  wall  pressures  versus  time  from  four 
calculations  for  9-kg  charges  in  a  1.83-m  radius  cavity. 
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Fig.  2  -  Conrparlson  of  wall  pressures  versus  time  from 
calculations  for  90“!^  charges  of  Pelletol  in 
a  4.57"ni  radius  cavity. 
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Equation  of  state  of  Pelletol  tased  on  Jones  and  Miller 
adiabatic  data. 
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Fig.  6  -  Wall  pressure  and  vail  pressure  average 

versus  time  for  a  450*'kg  charge  of  Pelletol 
detonated  in  a  1.83-m  radius  evacuated  cavity. 
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Fig.  10  -  Pressure  (bars)  versus  radius  at  indicated  time 
during  explosion  of  a  9-kg  charge  of  Pelletol 
detonated  in  a  1.83-in  radius  cavity. 


R  (m) 

Fig.  11  -  Density  versus  radius  at  Indicated 
times  during  the  explosion  of  a  9-1 
charge  of  Pelletol  detonated  in  a 
1.83-m  radius  cavity. 
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Fig.  17  -  Measured  and  calculated  radial  particle 
velocities  in  salt  for  a  250-kg  Pel let ol 
charge  in  a  4.57-®  raxiius  cavity. 
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Fig.  19  -  Particle  displacements  and  velocities  in 
salt  at  an  eighty  foot  (24.4-m)  range  for 
450-kg  Pelletol  charges  in  a  4.57-ni  radius  cavity. 
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experimental  and  theoretical  particle 
velocities. 


21  -  Power-lav  correlation  between  radiated 
and  explosive  energies. 


-67- 


Scoled  peak  Mach  number,  u^y/c/3  (per  meter) 
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Fig.  23  -  Theoretical  peak  particle  velocities  versus 


scaled  range. 
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Fig.  25  -  Theoretical  peak  radial  stresses  versus  scaled  range 
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Fig.  26  -■  Correlation  for  the  ex^jerimental  data  of 
peak  particle  displacement  in  the  salt  at 
various  distances  from  the  center  of  the 
cavities . 
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